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ABSTRACT 


This study, assuming stationary GausSian turbulence 
model, investigates the effect on the crossirg frequency of 
different spectral functions: von Karman, Kaimal and 
Teunissen for the x and z directions with five filters: 
ideal band-pass, ideal low-pass, "quadratic-type", "sine- 
type" andthe Hanning. The filters have a much greater 
effect than the spectral functicns. The estimated crossing 
frequency variation is as much as 50 percent among the 
guadratic-type, low-pass and sine-type filters. The Hanning 
filtering predicts crossing, rates up to thirty times, and 
the ideal one decade wide band-pass filtering predicts rates 
of eight times, higher than the ideal low-pass filtering. 
The variation, between the x and z direction, is less than 
10 percent for the von Karman spectrum, and over 40 percent 
for the Teunissen one. 
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I. INTRODUCTION 

Most flows occurring in nature and in engineering appli- 
cations are turbulent. In the earth's boundary layer, 
differential heating of the atmosphere produces pressure 
gradients, which are subsequently modified by the rotation 
of the earth, causing a complex velocity field. In this 
boundary layer the wind speed decreases as the surface is 
approached due to both the frictional drag of the surface 
and the drag of all bodies protruding into the air flow. 
These retarding forces are transmitted through the layer by 
shear forces and by the exchange of momentum due to the 
vertical movement of the air. The process of momentun 
exchange between layers is the mechanism leads to the gener- 
ation and decay of eddies which are termed turbulence. The 
resulting mixing of the air produces, along all three 
orthogonal axes, fluctuations in wind speed, commonly called 
gusts, which vary in size in both time and space. 

The effects of this atmospheric turbulence has’ been of 
continuing concern tothe aircraft or structure deSigner. 
Typical turbulence related prcblems are: the effects of 
turbulence on the fatigue life of the structure; the 
performance of control systems in turbulence and the deter- 
Mination of ultimate structural strength reyguired to accept 
unsteady loads induced by turbulence. In an attempt to solve 
these problens, a number of statistical models of turbu- 
lence, which endeavour to describe the turbulence in terms 
of as few parameters as possible, have been proposed since 
early in this century. To analyze these problems, one needs 
to know the crossing frequency, which is the frequency with 
which arandom function( here the filtered wind speed) 
crosses a prescribed value defined by a certain averaging 


time, say the hourly mean. 


The crossing frequency depends largely on whether or not 
the random function is filtered, the type of filter and toa 
lesser extent,the spectral density function. fhe effects of 
type of filter used and spectral density expression will be 
discussed in detail later. The spectral density functions 
have been proposed by Dryden [Ref. 1], von Karman [Ref. 2], 
Kaimal [Ref. 3] and Teunissen [Ref. 4]. The Dryden expres- 
Sion makes mathematical analysis simple but is much less 
accurate than the others. The von Karman spectral function 
was recommended by the Engineering Sciences Data JUnit 
(io «)Deulie) (Ref. 5] and this expressions was originally 
postulated for isotropic turbulence. The effect of departure 
from isotropic turbulence near the ground is allowed for by 
the variation with height and surface roughness of the 
appropriate variance 6 and length scale parameter L; since 
they typify the intensity and size of eddies constituting 
turbulence. The Kaimal expression, which is freguently 
used, was obtained recently for the surface layer over flat, 
uniforn, relatively featureless terrain in Kansas. The 
Teunissen model, which 1s a modified Kaimal one, was 
obtained for the generally rougher gross features of the 
upstream terrain. 

in the following discussion, we will review turbulence 
models and spectral density functions. The spectral density 
S(n), in non-dimensional form, is a function of the dimen- 
Slonless parameter nlj/U, where n is frequency and Il; is the 
length scale of the turbulence and U is the mean houriy wind 
speed. E.S.D.U. { Ref. 5] gives empirical expressions for 
the length scale obtained by analysing its coliection of 
world-wide turbulence data. Finally, for given values of 
nig/U, we will estimate the crossing frequencies with von 


Karman, Kaimal and Teunissen spectral expressions. 
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II. TURBULENCE MODELS 


aS Se Ee ae 


A. DISCRETE GUSTS MODEL 


Powell and Connell [Ref. 6] defined gusts as constitu- 
ting any series of discrete velccity-time events that can be 
defined froma turbulence time series according to some 
extrinsic criterion. There are two fundamentally different 
treatments of gust time; one is that the gust time is arbi- 
trarily fixed and the other is that the gust time varies. 
Examples of this latter type of gust definition, shown in 
Figure 2.1, represent a sample of wind fluctuation in the 
form of a wind-component time series. 

The GUSTO and GUST1 models define gust events in terms 
of gust amplitudes and characteristic times. Definitions of 
ampiitudes and times differ for the two models. The GUSTO 
model can be completely specified in terms of a poSitive 
peak amplitude +A,, and the time interval +T, between zero 
crossings on either side. It can also be specified in terms 
of a negative peak amplitude -A, and the corresponding time 
interval -T . Both definitions are expected to yield 
Similiar statistics. Powell and Connell combine both posi- 
tive and negative values into one set for their analysis. 

The GUST1 model can be specified in terms of a positive 
amplitude +A, which is the peak-to-peak amplitude between 
adjacent minima and maxima, anda time interval +T,; between 
the minima and the maxina; the positive sign indicates a 
positive rate of change in the variable. A comparable defi- 
nition can be made in terms of amplitude -A, and time -T,, 
where the negative sign indicates a negative rate of change 


in the variable. 
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GUST 0 Definition 


GUST 1 Definition 





Figure 2.1 Definitions of GUSTO and GUST1. 


The discrete model treats turbulence as a series of 
isolated gusts and the basic assumptions are 
1) Atmospheric turbulence can be modeled as a collection 


of isolated gusts randomly distributed along the 
structure. 


2) Gusts have random magnitude but fixed shape. 


3) The structure is a deterministic linear system with 
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sufficient damping, that each encounter with one of 
these gusts results ina single significant response 
peaks. 

This model is typically used to estimate ultimate 
strength requirements, however, it 1S not a satisfactory one 
for all aspects of the design problem. For example, its 
application in areasS such as structural fatigue, control 
system performance, or even extreme responses which involve 
lightly damped modes is highly questionable. In most 
instances turbulence does not occur as discrete gusts but as 
a continuous random disturbance. Furthermore, this model 
neglects most of the dynamic characteristics of the vehicle. 
Thus, although the necessity of evaluating the responses of 
a proposed structure to discrete gusts is still recognized 
aS an important part of the design procedure, the use of the 
discrete model to calculate most response statistics has 
largely given way in recent years to the use of turbuience 
models which attempt to take both the dynamic characteris- 
tics of the structure and the continuous nature of turbu- 
lence into consideration. This is known as the continuous 


Or power spectral model. 


B. CONTINUOUS MODEL 


The tasic idea of continuous model iS that atmospheric 
turbulence can be represented by a continuous stochastic 
process which acts asa disturbing influence on the struc- 


ture. The principal assumptions for this model are 


1j Each encounter of an structure with continuous atnos- 
pheric turbulence can be modeled as a deterministic 
linear system perturbed by three independent stati- 
Onary stochastic processes, which represent the long- 
itudinal, lateral and vertical gusts components. 


2) The spectral density of each random process has a 
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form characterized by three parameters, the gust vel- 
ocity standard deviation (6), the scale length(L), and 
the averaged speed (U). 

3) Each of the three gusts components is a Gaussian 

process. 

The assumptions of the Gaussian model make it possible 
to calculate tae statistics of any structure response for 
each region of turbulence as functions of the parameters L 
and 6 of that region. Thus, the assumptions of structure 
linearity and the Gaussian nature of turbulence permit their 
evaluation with a minimum of difficulty. These results will 
be dependent upon the assumed values of L and o along with 
the characteristics of the structure dynamics. 

The continuous model is used primarily to evaluate 
response statistics for selected conditions in continuous 
turbulence. For this application both the scale length and 
standard deviation of the turbulence as well as the struc- 
ture characteristics are fixed at values representative of 
certain conditions which could reasonably be expected to 
occur in service. 

The three assumptions made above need some reconsidera- 
tions. The first assumption stated that the turbulence is 
stationary and homogeneous and the dimensions of structures 
are much smaller than the scale lengths of the gusts compo- 
nents. These conditions are not always satisfied. The second 
assumption of the power spectral model is valid at low alti- 
tudes, but is not so certain at high altitudes, the problen 
made being more difficult by a lack of data at the latter 
altitudes. The third assumpticn stated that the turbulence 
is a Gaussian process. However, some experimental data indi- 
cate that the Gaussian turbulence model significantly under- 
estimates the frequency of occurence of high gust velocities 
[Ref. 7]. Furthermore, peak gust velocity cumulative prob- 


abilities of exceedance observed in atmospheric turbulence 
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behave iike exp(-x) rather than exp(-x*), aS predicted by 
the GauSSian distribution [Ref. 8]. Recently Reeves, et al. 
[Ref. 9] have discussed in great detail a non-Gaussian 
model. On the other hand, References(6) and (10) indicate 
that filtered wind speeds(particularly with a band-pass 5-50 
Hz filter) is reasonably Gaussian in character. In the 
present work, the Gaussian model is used exclusively. 

The statistical guantities of most frequent interest are 
the mean, variance, probability density, autocorrelation 
function, power spectral density and crossing freguency. The 
power spectral shapes usually assumed are those proposed by 
von Karman, Kaimal and Teunissen. The crossing frequency 
with these three spectral expressions will be studied in 


more detail later. 


{es 


III. SPECTRAL DENSITY AND PROBABIIITY DENSITY 


A. SPECTRAL DENSITY 


In the field of mechanical and electrical osciilations, 
a large number of phenomena are governed by equations of the 


forn 


2 
n oF + c St + ky = £(t) 


where, m, candk are constants; f(t) 1s a given forcing 
function and y(t) represents the response of the system as a 
function of the time t. 

If the above time dependent force f(t) is periodic, it 
can be expressed in terms of a Fourier series; this series 
represents the sum of a large number of forces of different 
amplitudes and frequencies. All elastic bodies have natural 
frequencies, and if this force is applied to an elastic 
system, andthe amplitudes of the components near a natural 
frequency are not very Small, then this force will drive the 
body into resonance. The body will then go through a very 
large number of stress cycles ina short period of tine, 
possibly leading to fatigue failure. 

To have stationary properties, a random Signal must be 
assumed to continue over an infinite time, and in sucha 
case neither the real nor the imaginary part of the Fourier 
transform converges to a steady value. In this case, we use 
the spectral density function, which has no convergence 
difficulty and which is applicable to a whole class of 


Similarly generated functions. 
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Spectral functions of atmospheric turbulence provide 
information on the freguency distribution of the kinetic 
energy of the various fluctuating velocity components. Used 
in conjunction with certain transfer functions, they provide 
information about the dynamic loading on, and response of, 
building and aircraft structures in the atmospheric wind. 

A considerable number of measurements of the power 
Spectra of gust velocities in the near-neutral atmosphere, 
or strong wind conditions, at varying heights and for 
different terrains are available. The von Karman spectral 
expressions, which were exclusively ised by) E-S-0-U - 
[Ref. 5], for each velocity component(x,y and z direction) 


are given by 


pnSu(n) _ Wnu (3.1) 
2 i: ~2 5/6 

ou (1+70.8i4) 
nSy(n) _ _4fy(1+755. 2fv) (eo) 
o2 (1+283.2n2)% ; 
nSw(n) Siw (1+755. 2hw) (B23) 
or (1+ 283. 2%3)'/ 

where Ny,= Lyn/U, h,= Lyn/U, hR,y= L,n/U and Ly; is integral 


length scale which is a function only of the terrain rough- 
ness and height above ground for the x and y directions; for 
the z direction it is a function of elevation only. oF is 
the variance of the wind velocity fluctuations about the 
hourly mean. These equations are in common use in the engi- 
neering literature and have the advantage that the integral 
Peoijeiescatoeei, sare treated as ‘iree* scaling parameters 
which are chosen to match the estimated scales for a partic- 
ular height and terrain type, while maintaining constant 


spectral shape. 
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Another commonly uSed Spectral expression is the Kaimal 
one; [Ref. 3] 


nSz;(n) __ 0.164 (f/f-) (3.4) 
o;? 440.164 (£/£, J” 
where nS(n) = logarithmic power spectral density 
Of = variance of i 
if = reduced frequency 
if = reduced frequency at the intercept of the 


extrapolated inertial subrange slope with the nS; (n)/ 
Oo; =1 line; or in our notation, foj =0.041Z/1; . 

The Kaimal expression is mere common in the meteorolo- 
gical literature and was obtained as a best-fit to surface- 
layer measurements over uniforn, flat, relatively 
featureless terrain. 

Teunissen [{Ref. 4] observed that von Karman model 
appears to te somewhat better than the Kaimal model at low 
freguencies. The Kaimal model largely underestimates the 
spectral content at these freguencies for all compo- 
hents(x,y,and 2). On the other hand, the Kaimal model shape 
appears to be better than that cf the von Karman nodel. 

Teunissen has proposed modified Kaimal spectral expres- 
sions, given by 


nSu(n) _-___105¢ (3.5) 
u,” (0.44+433£) 
pSy(n) _ 17f (3.6) 
> 
a? (0.38+9.5£) 
NSw (Nn) = 2E (3.7) 
U 2 0.444+5,3£% 
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where f=nL;/U and uy, 1S friction velocity. These expres- 
Sions were derived from the data which were obtained from a 
horizontal array of tower-mounted propeller anenone- 
ters (z=11m) during a five-hour period for which the mean 
wind direction was virtually perpendicular to the main span 
of the array, in the city of Toronto, Canada. 

As for the reason for the-higher spectral energies found 
by Teunissen in the low-frequency region of the spectra in 
comparison with the Kaimal valueS which was measured in 
Kansas, this is not entirely obvious, apart from the rougher 
terrain for Teunissen experinent. Teunissen has suggested 
that the power spectra over a wide range of terrain tyres 
May possitly be represented by uSing above expressions with 


appropriate values of such ‘terrain scaling parameters. 


B. PROBABILITY DENSITY 


In general, a knowledge of the spectral density ofa 
random process does not enable us to determine its prob- 
ability density or distribution. However, a common assump- 
tion which is reasonable in many applications is that 


atmospheric turbulence 1S a ‘normal or GauSsian process, as 


discussed above, with a probability density function for 
which 

: 1 ae 2 

(j= CX P| =1. 72673 Sis 
where i=u(t), v(t) or w(t) gare wind-speed fluctuations, 


which have a zero mean. Note that Equation (3.8) depends only 
on the standard deviation of the gust velocity. Thus, if the 


probability distribution of o0 is known, the distribution 
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function of the gust velocity can be found by integration, 


aos ae 


<p 
P - frviyas be) 


In practice, atmospheric turbulence contains patches of 
a Significantly non-Gaussian nhature (particularly in the 
lower 30m) when larger gusts and longer lulls may occur more 
frequently then indicated by the Gaussian distribution. 
Kaimal et al. {Ref. 10} observed that the assumption of a 
Stationary Gaussian process has validity only for the 


filtered data, not for the unfiltered time series. 
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The instantaneous horizontal wind speed is considered to 
be composed of a mean longitudinal component J plus’ fluctu- 
ating component vw;, where the subscript i designates the x, 
Ye Oxr Zz Component ina coordinate system with x oriented 
along the mean wind vector and Zz upwards. Then the wind 


speed as a function of time is given by 
V(t) = U + v;(t) (4.1) 


The average frequency of fpositive-slope level crossing 
Noi, with which the wind speed fluctuations exceed their 
zero value, for a GausSian distribution of turbulence, is 


given by 
1 / 


Ne2 Ne2 
No; -| [otrarscteran J frims; (n) a (4.2) 
n 


Ch Net 


where n,. is cut-off frequency, iis u,v,or wand F(n) isa 
filter function. 

The value of No; depends on the mean wind speed U, L; and 
turbulence intensity o7/U, which is a measure of the magni- 
tude of turbulence fluctuations and defined as the ratio of 
the standard deviation of the instantaneous fluctuating 


velocity component to the mean wind speed averaged over 


arbitrary time, usually hourly, and depends only on the 
height 2Z and the terrain roughness Z, , in a_ stable 
atmosphere. 
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A. FILTER FUNCTIONS 


It is clear that the guality of the crossing frequency 
estimates obtained by means of Equation (4.2) will depend 
strongly on the form of the filter function F(n) and, hence, 
on the selection of the time function h(t). Thus we will 


consider some examples of the function h(t) and filter func- 


tion F(n) corresponding to theug. Let the time function be 
given by 
h, (t) = 2 Sing2e7 7) (4.3) 
T Zie7 L 


If we take Fourier transform, the Equation (4.3) becomes 


Sn Be eae 1 when |nj{ < 1/T (4.4) 
1/2 when jn| =1/T 
0 when j[n{ >1/T. 


Graphs of the functions h,(t) and F, (n) are shown in Figure 
4.1 and Figure 4.2 respectively and they represent an ideal 
band-pass filter. 

The simple impulse response function h,(t) is the ordi- 


hary moving average: 


h, (t) {00 when — T/ 2 Sue 27 2 (4.5) 
0 when Jt] > T/2 


The corresponding filter functicn F, (n) (F, (n) =H, (n)*), where 
H,(n) is the Fourier transform of the function h(t), has 
the form 


2 
F,(n) = ( sintTn (4.6) 
iin 
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Figure 4.1 Time Function h, (t), hy(t), hs; (t) and h, (t) - 


The function h, (t) and F,(n) are shown in Figure 4.1 and 
Figure 4.2 respectively. The figure shows that F,(n) is a 
low-freguency filter with cut-off frequency ne =1/T. The 
Shape of this filter shows that frequencies above n are 
transmitted nonuniformly and those below n. are completely 
suppressed. : 

The impulse response function h, (t), called the smoothed 
moving average is: 


h,(t) =(1+cos(2™/T)t when -T/2 <t¢ T/2 (i a7 
it 
0 when eines 1/72 
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FILTER FUNGHON 
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Figure 4.2 Filter Function F, (1). F2(2), P,(n) and F, (n)- 


The corresponding transfer function F, (n) has the forn 


F, (mn) = ( Sines) nee (4.8) 
Mie | i. 
Tn (n*- (rtn}| 
The advantage of F, (0) over F, (2), is that it has much 


Smaller side maxima. 


Another commonly used filter function is the guadratic 
filter function given by 


Fy (n) = 1 (4.9) 
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One can find the impulse function of Equation (4.9) by 
inverse Fourier transforn, i.e., 


o 


(ng +n? ) 
- 


1 


Ya j2unt 
IT 2 dn 


hy (ty 


st ((n2 +n? M7 *ccs2nntdn + 


ene 2 +n? )? 
Sa jing *n (jsin2mnt) dn (4. 10) 


The second term of above equaticn is zero, since, it is an 
odd function. Thus Equation(4.10) becomes 


o 


“V4 
h, (t) =1- five +n*)* cos2nnutdn 
© 


or hy (t) = Ko (nt) (om) 


where K,is the modified Bessel function. 

The functions h, (t) and Fy, (n) are Shown in Figure 4.1 
and Figure 4.2. 

The cut-off frequencies are defined formally as_ the 
frequencies at which the filters reach the first zero. But, 
when a real filter is used, we must find some equivalent 
cut-off frequency, which is defined as the cut-off frequency 
of an eguivalent ideal filter whose area is equal to that of 
the real filter. [Ref. 11] Greenway [Ref. 12] and Powell 
and Connell [Ref. 6] have used the value at the half- 
power (3db) point of the filter as the cut-off freyuency. 
Vinnichenko et al. [Ref. 11] defined the equivalent cut-off 
frequency as 


4) 


Ne = [rena ; (4. 12) 
Oo 
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Figure 4.3 Effects of Cut-off Frequencies. 


According to the Equation(4.12), the cut-off frequencies 
depend on the filtering functicn. The effects of cut-off 
frequency on the crossing frequencies are shown in Figure 
4.3, which shows that the higher the cut-off frequency leads 


to lower zero crossing rates. 
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B. CALCULATION OF CROSSING FPRECUENCIES 


In this section the frequency of zero crossings are 
estimated using von Karman, Kaimal and Teunissen spectral 
density functions with various filter functions which we 


have discussed above. 


1. Using the von Karman Spectral Density Function 


Introducing the von Karman Spectral density function 


for the x-direction S,(n) into Fquation(4.2) yields 


I 
Yica n 1: 
Nox = fz F (n)dn / “F (n) dn = (4.13) 
(1470.8n2)/ ) (1470. 81254 
Nc Ne; 
For numerical calculation, set y=n/n,g, and D=n,., /10N¢, and 


Equation (4.13) transformed to 


10D 


ly, 

BY. ie 2, 

Nox =Ne, y* dy i dy (4. 14) 
(1470 8g, y7y4 (1470. 8ii<, y?7)@ 


where Nae is the dimensionless frequency of positive-Sslope 
level crossings for the x direction Noxly/U. Equation (4. 14) 
is for the ideal band-pass filter function of D decades 
filter width. 

The upper integral of Eq. (4.13) will not converge as 
Neg goes to infinity because the empirical expression for 
the turbulence spectrum employed does not account for the 
high-frequency range where viscouS dissipation damps out 
turbulence fluctuations. One can achieve convergence by use 
of a suitable filter function. Introducing Equation (4.9) 
into Equation(4.13) yields 


0 y, 


J (1#70.88Z y7)% (1+ y?) (1470. 872 y? {1+ y?) 


oO 
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If we take Fr, (n), here called the sine filter function, 
instead of quadratic filter, with the cut-off frequency 
n.=0.44/T which Greenway {[Ref. 12} has used, and substitute 


into Equation(4.6), we find 
2 
EA (n)=sin*(0.44qn/n,)/(0-44m D/n,) - (4. 16) 


Introducing Equation (4.16) into Eguation (4.13) and setting 
y=0.44nn/n, yields 


m iy, 
sin* (y) dy 
thy le (1+70.8(He /0-.44m)? y?| (4.17) 


0. 44 fe (sin? (y) /y?] dy 5 
[1+70.8 (i /0. 44m)? y?| : 


a O 


If we use Ps (n), which is the Hanning frequency 
window, and equate cut-off freguency nh. =0- 3957 Peas derivad 
by Vinnichenko et al. [Ref. 11], then substitution into 


Equation (4.8) gives 


2 
F, (n)= ( sin (0.3751n) 1* (4. 18) 
0.375nn /[m?-(0.375nny] 


Substituting Equation(4.18) intc Equation(4.13) and setting 
y=0.375an/ne and =, /0.375N yields 


| 
z fs sin* (y) dy ie 
Nox =%, 4 A(14+70.8%o y2)% (n?-y2)? (4. 19) 


(sin? (y) /y?]} dy 
J(1+70.8 42 y2y4 (W-y?)? | 
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TABLE I 
Crossing Frequencies(von Karman x-direction) 


t 


Te IDBP QUAD SINE HANNING 
0.100 0.3972 Gactas 0.0865 17-4480 
Oeg125 0.4838 0.1116 0.0955 21.8074 
0. 150 025709 O12 5.0 0.1117 26. 1670 
0.175 0.6585 Os 384 OO, 1232 Soe oO 2on 
0.220 0.8172 0.1604 0. 1424 38.3741 
0.500 1.8193 O2-2737 0.2391 87.2049 
1.000 3. 6255 0.4298 0.3699 174.4062 
2.000 Teoagh 0.6765 0.5744 348.8093 
5.000 18.0931 e236 5 10383 872.0234 
10.000 36.1847 [ges bee pe. 1.6267 1744.0481 
15.000 54.2766 265565 221230 2616.0669 
22-000 79.6055 3.2989 Zetl 30k 3836.9021 
50.000 180.9212 5.6928 4.6950 8720.2461 
75.000 271.38 16 Tat 5.5 5 6.1414 --- 
100.000 361.8425 9.0288 7.4326 --- 
150.000 542.7637 11.82€8 9.7292 --- 
175.000 633.2241 13.1054 10.7787 --- 
200.000 723.6853 14.3243 11.7794 --- 
225-000 814.1465 171.4933 12-939 2 nee 
250.000 904.6067 16.6196 13.6641 --- 


he =cut-off fre uency 
IDBP=ideal ban = Fass fit Ler 
QUAD="quadratic-— YP e" filter 
SINE="sine- ane iiter 
HANNING=HANNING filter 
IDLP=ideal low-pass filter 


haa 


| 


ne nn en ee ee | 


The results of numerical calculations of Eqs.(4.14), (4.15), 
(4.17) and (4.19) are given in table I. The plots for above 
equations are shown in Figure 4.4. The plots showing how 
well these crossing frequencies co-incide depend on the 
length scale, mean wind speed and filter functions. 

The von Karman spectral density function for the 
z-direction is given by Eguation(3.3). The equation of the 
crossing frequency for the ideal band-pass filter function 
can be obtained by SUbSertucang Hquation (3. 3) into 


Equation (4.2) and is given by 


Zo 
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Figure 4.4 Crossing Frequencies with von Karman 


Spectral Expression for the x-direction. 
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Using the same procedure as for the x-direction, the 
z-direction crossing freguency for the quadratic filter 


function is found to be 


1; 
oo ae Z 
il (1+755. 2, y?) y2 dy 
Neea=n, (14+283. 2% a cea? (4.21) 
(14755. 22 y?) dy 
| (14283. 2n2y? )4(1+y2) 


and, for the sine filter” Eunction, i: one defines 
B./0.44 =", the crossing fregquency for the z-direction 
becones 


ira) 1 
755. 2 (Ne a= (y) dy 
=, Ree283.2%.ya-° 4 (4. 22) 


fi +755. 2 (Ney) J (sin? (y) /yJay 


(1+283.2%2y2)'% 


Introducing filter function Fz (n) into Equation(4.20) and 
setting y=0.375nn/n, and fic /0. 375N=%, yields 


oa Y/, 
(1+755.2 (%y)] sin? (y) dy 


Nox Me 4 ~2_(1#283. 22y?)S (me -y?)" (4. 23) 
(1+755.2(%e Jt sin? (y) /y’Jdy 
(14283. 24 2y2)% m7-y?)? 
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ee 


TABLE il 


Crossing Freguencies(von Karman z-direction) 


: 
| 
| 

De IDBP QUAD SINE HANNING 
| 0. 100 0. 3342 Opal 0.0816 17.4448 | 
| O- 125 0.46 86 0-186 0.0934 21.8044 i 
0. 150 0.5567 0. 150 0.1044 26.1641 | 
0.175 0.6454 0.163 0.1144 50 .9 Zoic i 
0.220 0.8059 Os 185 0.1320 SO- 37 16 | 
0. 500 1. Sidley 0.297 Oe218e 87.2036 
1.000 3.6207 0.452 O=5575 174.4053 
2-000 71.2380 0.696 0252338 348. 8093 
5.000 18.0925 1.248 Dig slit sye) 872-0229 | 
10.000 36.1844 1.951 1.4870 1744.0461 
15.000 54.2764 22530 1.9407 2616.0693 
22-000 79.6054 3.246 22-4978 3836.8989 
50.000 180.9214 5.493 4.2994 8720.2266 | 
715.000 2/71. 3823 Fig Ay 5.6256 --- | 
100.000 361. 8430 8.627 6.8094 —— 
150.000 542.7634 11.129 8.9150 ——— 
175.0009 633-2250) 2.53 0 9.8773 --- 
200.000 723.6838 13.474 10.7948 -<-- 
225.000 814.1440 14.572 11.6748 a 
250.000 904.6057 15263 C lZeo228 --- 

feos ae 2 


The results of numerical evaluations of Eqs. (4.20), 
(4.21) , (4.22) and (4.23) are listed in the Table II. The 
plots for above equations are shown in Figure 4.5. 

For comparison, the crossing frequencies with von 
Karman spectral expression for the x and z directions are 
plotted in Figure 4.6. The crossing frequency for the z 
direction is less than the x direction except ideal band- 
pass filter case, in which the frequencies are identical for 
both directions. 

The crossing freguency plots presented in ‘Figures 
4.4 - 4.6 offer a detailed view of how the crossing 
frequency depends on the filter function. Thus, in designing 
a system, the type of the filter has to be chosen carefully 
to match the behavior of the particular systen. 
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Figure 4.5 Crossing Frequencies with von Karman 
Spectral Expression for the z-direction. 


2. Using the Kaimal Spectral Density Function 


From Equation(3.4), the asymptotic spectral behavior 


can be written as 


: “3 
nSz(n) _ (f/f, ) : £>>f, (4.24a) 
o;* 0.164(f/£,) 3; £<<£, (4.24b) 
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Figure 4.6 Crossing Frequencies with von Karman 
pectral Expression for the x and z-direction. 


The expression for the low frequency behavior of the 
spectrum is given 


Set) = ( 4 Li/z)F 
oF (4. 25) 
where L; represents the length scale. Eguating 


Equation (4.25) and Equation (4.24b) becomes 
Le; = 0.041Z/ (£,)x (4.26) 
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HOE sthne x-direction only 


Lu = 0.041Z/fox 
Or foy= 0.041Z/Ly (4.27) 
Substituting Egquation(4.27) into (3.4) yields 
Nowe. 0.164 (£/0. O41Z/Ly) 
cite 140.164 (£/0. 04 1Z/Ly) 
or nSy (n) un 
o., 1+33.641n 3 (4. 28) 
where. Tl hhLy/U. Then the crossing frequency for the 
X-direction becomes 
72. 
as 0 2 ®” 
Nox = F(n)n*dn i} F(n) dn 
}(1+33.6410 4) }( 1+ 33.6418) (4.29) 
Substituting the quadratic filter function into 
Eguation (4.29) yields 
y, 
o 
[ n* dn 2 
Nox = |~o(1+33. 641874) (ng +n?) (4. 30) 
O° 
rf dn 
J(1+33.6410%) (ng +n?) 
Replacing n/n, by y in Equation (4.30) leads to 
1 
b /; 
mw) y7dy 2 
en ~ oe PY 
Nox =D, s[1+ 33.641 (fe y)73] (1+ y2) 
dy (4. 31) 


6 


lis 


641 (fic y)?] (1+y?) 
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Setting y=e" in Eguation(4.31) gives 


o /2 
: ei% ax 
Nox =Ue 2 4% [1433-641 (fc e*yA] (142%) (4. 32) 
id e* dx 
) (1433-641 (7 eX 73] (1+e2%) 
Pt we introduce the sine filter funet ion into 


Equation (4.29), the crossing frequency becomes 


Vo 


- air > 
f sin* (0. 44nn/nc ) n? dn 
N = 0 


(14+33.6418 7) (0.4ut/n,)* n* (4. 33) 
ia Sina 024 on, ca an 
) (14+33.641n7°3) (0.44m/n, )? n? 


ox 


Substituting 0.44mn/n- for y in Eguation(4.33), and define 
fey OaiW t= remy wedldic 


fi sin? (y) dy 
=. | % 1433.641 (ty)? (4. 34) 
“{sin? (y) /y?] dy 
f 14+33.641(7, vy) 


Ze 


OX 


For the Hanning filter, substituting Equation(4.18) into 
Equation (4.28), setting y=0.375,tn/ne and define 
i, /0.3751=Ra, yledds 


36 


2 
| (4.35) 


Crossing Freguencies (Kaimal x-direction) 
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The results of numerical calculations of Eqs. (4.29), 


in Table III and the 


are given 


(4.34) and (4.35) 


(4.32), 


plots for above eguations are shown in Figure 4.7. 


that, 
the difference in 


and Eq. (4.26) 


Ege (3-4) 


for the Kaimal spectral density function, 


It is clear from the 


and z directions is only 


crossing freguencies retween the x 


Thus the plots for the z-direction 


in the length scale "L;. 


are exactly same as the x-direction. 
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Figure 4.7 Crossing Frequencies with Kaimal 
Spectral Expression far the x-direction. 


3. Using the Teunissen Spectral Density Function 


The modified Kaimal spectral density expressions, 
which were derived by Teunissen, are given in Eqs. (3.5), 
(3.6) and (3.7). These eguations were derived by adding a 
scaling parameter to each of the Kaimal spectral equations 
and evaluating these parameters from the observed data. 


These equations are function of friction velocity Uy - In 
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order to estimate the crossing fregyuency as we have done 
earlier, we need to change that eguations as a function of 
waietances) In £.5.).U. (Ref. 5], the relations between 6; 


and uy are defined as 
0.460; /u,=(0; /U) log (Z/2z,) (4. 36) 


Teunissen derived the average values of o;/u,, given by 


oy ee Day (4.37a) 
a7 ly el 2 (4.37b) 
O,y/Uy, = 1.27 (4.37c) 


From Equation(4.36) and Equation(4.37), one can get 


Us (u-component) = 0.124 on (4.3 8a) 
uy (v-component) = 0.271 oy; (4.3 8b) 
Uy (w-component) = 0.620 oe (4.38c) 


Siaesilostrtuting Eq. (4.38) into Eys. (3.5), (35,60) 9.and “3.7) 
yields 


nSyu(n) | 13n (4. 39) 
on? = ~5/s 
- (0. 444+33n) 
nSy(n) 4.61n (4.40) 
o,,? (0. 38+9.5n)4 
pSw(n) _  __1624n (4.41) 
a2 0.4445.31% 


Using methods similar to the akove, the crossing frequency 
for the x-direction, with the ideal band pass filter, 


becomes 
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rs) 
a n? dn 1s dn 
(0.444+33n)?/ J (0.44+33n) 
Net 
p | 


MY, 
ox = Neg loneage = a ee aa 
(0. ae Jo Nee Gize wi” 


If one uses the sine filter function, With cut-off 


ma 
1 


or 


frequency n.=0.44/T, the crossing frequency becones 


90 1 
| sin? (0.44nn/n¢) n? dn /2. 
N — 


6a (0. 44+33n) (0 .44mn/ne rg (4.43) 


|e Sin? (0.44nn/n.) dn 
J (0.44+330)"3 (0.44nn/n, )? 


Setting y=0.44nn/n, into Eq. (4.43), and define i, /Oe uate 
yields 


go WZ 
[ sin? (y) dy (sin? (y) /y2Jdy )* 
i ox a eee oo.  —— ae 
Nox =e ; (0. 44+337. yy? : (0.444+33% y)% (4. 44) 
For the quadratic filter function, Nox becomes 


yy 
y7 dy - 


Nox = De ee y)3 (1+y?) (4.45) 


ieaue 
J (0.44+338, y)¥% (1+y?) 
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TABLE IV 
Crossing Frequencies (Teunissen x-direction) 


— 


UW 
GI 





FADNOCANAMTCADOKElHNOCCMWMA 
ONNTFOKnTMOKOTDONMNE™Nnaowr 
MON FID OW DONT ODEeMaAMow 
DOCCOCEKNFROMAHOMAM NWO 
eer ee @ @ @ @ © @ @mUcheemUCc HhUCUcOFrmhUhCUC OhUCUMHhUCUc HMhCUhhUCUCFhlhUh 


CODTOOCOCOCOTEKNANAM=ANWMWiw 


NFHS FREMNNAMNODANEHOTONe 
NAINA OASWIDOIFNCCOrPhIMKMNe 
PFPIFEOM ANTM CONCTDOODONDOOCO 
DODCVCOTTrMNATNOMONAONAM 
ee 8 @s#8s 8 @¢@8e¢ « © © 6.8hUmFmUCUC HMmLUC HhUrhhUCUrPhUlUrlhUlUrhUlCUMD 


COCOOCCOOCOTe NM ANWMWwWOr 


Oo 
BOYD ONDNDO ODOM TAI NOWM Stine 
WOM ODWWONWTOMNRDMOWMONN 
CODOOSTNAMNMNOMWOMNNOOL) -« 
eeevsee5r2:etsee Gees © «€ © @« & 0 


OODVTOVCOCVCOFTENMNSTUDOAE” 


TmmOnToOWMO 

MAMAN MMAFAWNWworvww 

NO FHDOTATOMYUIMOMAWUIMNAWN SE” 
NAMM OOCDPODOKHAMOMNAW rw 
MWOMNTONNATTRNO « © ee we ee 
MTWOOCBDMION «6 « 2 OrrnNnnnsast 
© ee © © 8 © wg POO THOM OANNANTO 
DODCOOTMM EMME KAMNOPON 


DOODVDOO°CoO 

COCDOTOCOCCOOO 
CMNOMNODOCOCOCOCOCOOCOCOCOCO 
ONION OOOTOOOO0C0 «€ © © « « «@ 
TM OMONINOOC «©« wo eo e OOMOMWMO 
ee © e 0 8 8 ke OMNOMOWMONW 
ODD OCOO CANN CCA Cerrar nN 


The results of numerical calculations of Eqs. (4.42), 


Table IV and the plots for 


are given in 
these equations are shown in Figure 4.8. 


and (4.45) 


(4.44) 


ex pression 


modified Kaimal spectral 


the 


for 
If we put ideal band pass 


The 
z-direction is given by Eq. (4.41). 


and spectral expression for the z-directon 


filter function 


into Eq. (4.2) it yields 
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Yic2 
dn 
0.44+5.3 


Net 


n 
0.4445.3n 7 
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or 


‘la 
Ne2 
kc, iP dy (4. 46) 
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Figure 4.8 Crossing Frequencies using Teunissen 
Spectral Expression for the x-direction. 


W2 


If one introduces the sine filter function, the equation for 


the crossing frequency becomes 


fe (2 <a 2 a 2 
Nog =| Yo(0.44+5.30%) (0.44un/ncy (4.47) 
ae (0.44nn/n, 
3 / 
[ Sin*~ (y)dy 
OL Nee = bh 0.444+5.3 (Ne yy? (4.48) 


“(sin? (y) /y?] dy 
, 0.444+5.3 (Ty)? 


If one introduces the quadratic filter function, the equa- 


tion for the crossing frequency becomes 


Y 
] 
w 
y* dy 
-1 A[0. 4445.3 (RL y)4] (1+y?) (4.49) 


= 
il 
oR 


a 
ifs 4445.3 (Ry)? J (1+y*) 


The results of numerical calculations of Eqs. (4.46), 


(4.48) and (4.49) are given in Table V. The plots for above 
equations are shown in Figure 4.9. 
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TABLE V 
Crossing Frequencies (Teunissen z-direction) 
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es 


the crossing frequencies 


For comparison, 


with the 


expressions for 


Kaimal and Teunissen spectral 


von Karman, 


the x-direction with various filter functions are plotted in 


crossing 


As shown in Figure 4.10, the 
Kaimal spectral expression are 


4.10. 
freyuencies predicted by the 


Figure 


higher than for the others. 


the crossing 


In Figure 4.11, 


and the 


spectral expressions 


The 


all three 


with 


frequencies, 


figure characteristics are 


are compared. 


Z-dl rection, 


x-direction, except that the 


the 


for 
crossing frequencies are lower. 


to those 


Similar 
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Figure 4.9 Crossing Freguencies with Teunissen 
Spectral Density Function for the z-direction. 


C. EXAMPLE CALCULATIONS 


In the previous section of tnis study, we derived equa- 
tions for the crossing frequency. The purpose of this 
section is to illustrate the application of these equations 
to a given turbulence condition, and to estimate the 
crossing frequencies for a low-pass filter. 
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Figure 4.10 Comparison of von Karman, Kaimal and Teunissen 
Spectral Density Function for the x-direction. 


The dimensionless crossing freguency, Noze is a function 
of D, (=nlyg /U). The length scale parameter Lg and surface 
roughness parameter Z, are given by E.S.D.U. [Ref.j 5]. For 
convenience these data are included here as Figure 4.12 for 
length scale parameters, and Figure 4.13 for the surface 


roughness parameters. 
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Figure 4.11 Comparison of von Karman, Kaimal and Teunissen 
Spectral Density Function for the z-direction. 


Suppose we have a site for which 


Z = 20m 
Ne= 50HZ 
U = 25m/sec and the turbulence condition is that 


for an off-sea wind in a coastal area. 
Wath this choice of turbulence condition, the surface 
roughness parameter z, can be found in Figure 4.13, and is 
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Z_e=0.001. The length scale parameters are given in Figure 
4.12, and are given as Ly=105, 1,=36 and Ly=6.5. 
The next step is calculation of the dimensionless 


cut-off frequencies, i.e., 


Toy= DLy/U = 210 (4.50) 
Dey= nly/U = 72 
Lega Blg/o = 1s 


The crossing freguencies are calculated by the use of 
the given crossing frequency vs cut-off frequency graphs, 
with dimensionless cut-off frequencies which are given in 
Ege (4.50). 

The final step is to convert the dimensionless crossing 
frequency toa dimensional freguency by multiplying the 


factor U/lj, dees. 


Nox] Noss uae (4.51) 
Noz = Noz U/lw 


' The results of calculation cf the crossing frequency for 


a given turbulence condition are given in Table VI. 
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TABLE VI 
Example Calculations of the Crossing Frequency 
von Karman (x-direction) 
REMARK oeEE footie QUAD SINE 
Nox 75.0 110'.2 14.6 12.0 
Nox 172057 2. 428 3.476 2.857 
fie 210.0, 0/Lu=0.238) 
von Karman (z-direction) 
REMARK ITBP IDIP QUAD SINE 
Noz 4.70 1.50 2.40 1.80 
Noz 18.076 5. 769 9.231 6.923 
fticg=13, U7 lw=3. 846) 
Kaimal(x and z-direction) 
REMARK IITBP IDLP QUAD SINE 
Nox 75.0 10.9 16.50 12.2 
Ox 17.857 2.595 3.928 2.905 
Oz 4.70 1.70 2605 205 
| Noz 18.076 6. 538 10. 192 7.885 
: Teunissen(x and z-direction) 
| REMARK IDBP oDLe QUAD SINE 
Nox 75.0 6. 08 9.30 6.50 
| Nox 17.857 - 448 Deo 1.548 
Noz 4.70 2510 B30 20 
Noz 18.076 8. 076 12.692 9.230 
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Figure 4.12 Values of Length Scale Parameter lL; 
(Reproduced from ESDU 74031). 
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Figure 4.13 Values of the Surface Bom gntess Parameter Z, 
(Reproduced from ESDU 74031). 
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V. CONCLUSIONS 

This study has estimated the crossing frequency based on 
the three spectral density functions, which assumes turbu- 
lence to be a homogeneous, stationary GauSSian process, with 
five different filter functions for the x and z directions. 

This estimation has shown, for the x direction, that the 
Kaimal spectral function gives about one percent higher 
crossing freguency than the von Karman, and forty percent 
higher than the Teunissen one, at Ne =100 for guadratic-type 
filtering; for the ideal low-pass and sine-type filters, the 
difference between the Kaimal and von Karman spectral func- 
tion is only a few percent, kut for the Teunissen one is 
More than forty percent than for the other spectrums. 

For the z direction, the Teunissen spectral function 
gives about ten percent higher crossing frequency than the 
Kaimal one and about twenty percent higher than the von 
Karman one, at n,=100 for the quadratic-type, sine-type and 
low-pass filtering. However, this difference becomes to zero 
aS non-dimensional cut-off frequency goes to zero. 

The crossing frequencies for band-pass and Hanning 
filtering are only a fraction of a percent different for the 
three spectral functions. 

The present study shows that Hanning and band-pass 
filtering predicts higher crossing frequencies than low-pass 
filtering. This difference is greatly amplified with the 
higher frequency. 

If some site statistics are available, as for the 
example calculation, then the number of crossings that the 
wind-speed fluctuations exceed their zero value per a 
certain period can be computed and an estimation of fatigue 


life can be obtained. 
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AS a conseguence of the results studied herein, itis 
suggested that, in future work, careful consideration should 
be given to the type of filtering used to model structural 
behavior and determine more accurately the fatigue life of a 


systen. 
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